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Abstract 



Nonlinear dynamics of a bouncing ball moving vertically in a gravi- 
tational field and colliding with a moving limiter is considered and the 
Poincare map, describing evolution from an impact to the next impact, is 
described. Displacement of the limiter is assumed as periodic, cubic func- 
tion of time. Due to simplicity of this function analytical computations 
^ ■ are possible. Several dynamical modes, such as fixed points, 2 - cycles and 

r*^ ' chaotic bands are studied analytically and numerically. It is shown that 

chaotic bands are created from fixed points after first period doubling in 
a corner-type bifurcation. Equation for the time of the next impact is 
solved exactly for the case of two subsequent impacts occurring in the 

CO . same period of limiter's motion making analysis of chattering possible. 

O 
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1 Introduction 

In the present paper we study dynamics of a small ball moving vertically in a 
gravitational field and impacting vifitli a periodically moving limiter (a table), 
j^ , This model belongs to the field of nonsmooth and nonlinear dynamical systems 

[U [31 ini S] . In such systems nonstandard bifurcations such as border-collisions 
and grazing impacts leading often to complex chaotic motions are typically 
present. It is important that nonsmooth systems have many applications in 
technology H [71 18] . 

In the bouncing ball dynamics it is usually difficult or even impossible to 
solve nonlinear equation for an instant of the next impact. We approached this 
problem assuming a special motion of the table. Recently, we have considered 
several models of motion of a material point in a gravitational field colliding 
with a limiter moving periodically with piecewise constant velocity [9l 110) and 
velocity depending linearly on time |11| . In the present work we study the model 
in which periodic displacement of the table is a cubic function of time, carrying 



out our project to approximate the sinusoidal motion of the table as exactly as 
possible but preserving possibility of analytical computations [12] . 

The paper is organized as follows. In Section 2 a one dimensional dynamics 
of a ball moving in a gravitational field and colliding with a table is reviewed 
and the corresponding Poincare map is constructed. A bifurcation diagram is 
computed for displacement of the table assumed as cubic and periodic function 
of time. In the next Section dynamical modes shown in the bifurcation diagram 
such as fixed points, 2 - cycles and chaotic bands as well as the case of N impacts 
in one interval of the limiter's motion are studied analytically and numerically. 
We summarize our results in Section 4. 



2 Bouncing ball: a simple motion of the table 

Let a ball moves vertically in a constant gravitational field and collides with a 
periodically moving table. We treat the ball as a material point and assume 
that the limiter's mass is so large that its motion is not affected at impacts. 
Dynamics of the ball from an impact to the next impact can be described by 
the following Poincare map in nondimensional form [T3] (see also Ref. [2] 
where analogous map was derived earlier and Ref. [15] for generalizations of the 
bouncing ball model): 



7ym+i) = 7y(T0-A2+i+A,+iV„ (2.1a) 



v,+i = -i?y, + 2i?A,+i + 7(i + i?)r(r,+i), (2.1b) 



where Ti denotes time of the i-th impact and Vi is the corresponding post- impact 
velocity while A^+i = Ti+i — Ti. The parameters 7, R are a nondimensional 
acceleration and the coefhcient of restitution, < R < 1 5 , respectively and 
the function Y (T) represents the limiter's motion. 

The table's motion has been usually assumed in form Ys{T) — sin(r), cf. 
[m [15] and references therein. In this case it is basically impossible to solve the 
Eq. (|2.1a|) for T^+i. Accordingly, we have decided to choose the limiter's periodic 
motion in a polynomial form to make analytical investigations of the dynamics 
possible. In our previous papers we have assumed displacement of the table as 
piecewise linear periodic function of time [SJ [TU] as well as quadratic [TT] . In 
this work we study dynamics for a cubic function of time Yf. (T): 



Yc (T) = 12\/3f (f - 5 j (r - 1 j , (2.2a) 

Yc (T) = 6\/3 (er^ - 6r + 1) , (2.2b) 



with T = T — [T\ , where [xj is the floor function - the largest integer less than 
or equal to x (i.e. < T < 1). 

Since the period of motion of the limiter is equal to one, the map (12. ip is 
invariant under the translation T^ — >■ T^ + 1. Accordingly, all impact times Ti 
can be reduced to the unit interval [0, 1]. The model consists thus of equations 
(|2.ip . (|2.2[) with control parameters R, 7. 
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Figure 1: Bifurcation diagram, R — 0.85, 7 G [0, 0.06]. 



In Fig. [T] above we show the bifurcation diagram with impact times versus 7 
computed for growing 7 and R — 0.85. It follows that dynamical system (|2.ip . 
(|2.2p has several attractors: two fixed points which after one period doubling 
give rise to chaotic bands and two other fixed points which go to chaos via period 
doubling scenario. There are also several small attractors. We shall investigate 
some of these attractors in the next Section combining analytical and numerical 
approach. General analytical conditions for birth of new modes of motion were 
given in jl6) . 



3 Analytical and numerical results 

3.1 Fixed points and their stability 

We shall first study periodic solutions with one impact per k periods. Such 
states have to fulfill the following conditions: 



14+1 = Vn = vi''/'\ T„+i =Tn + k = Ti''^^^ + k 



(A; = 1,2,...), (3.3) 



where: 



T, 



(k/ 



'^ e (0, 1) , v.j'/'^ > 7r,, (ri'=/^)) 



(3.4) 



Substituting these conditions into (j2.ip . (|2.2p we obtain two sets of fixed 
points: 



< T 



(fc/i) 



5-^V97^ + ^^fc7iTl<I 



(3.5) 



v: 



(fe/i) 



where the impact occurs in time interval TK^^ £ (O, ^) and 



I < riK)^ = l + ^V97^ + ^fc7i^<l (3-6) 

yi'-/') = k 

with impacts taking place in time interval TS^, e (i l). 

Solutions p.Sp fulfill physical requirements and are stable in the following 
interval of 7: 



fci-l < 7 < 54(TW (3^(^' - 1) + V9fc2(i?2 _ 1)2 + 12(i?2 + 1)2) (3.7) 



V3i 1-R 
18 '^ 1+fl 



where lower bound is a consequence of T^'/-, > while the upper bound follows 
from the condition that eigenvalues A of the stability matrix obey |A| < 1. Ac- 
cordingly, for R — 0.85 there are only four stable fixed points with k — 1, 2, 3, 4, 
shown in Fig. H]- they appear at T = and 7 = 0.0078, 0.0156, 0.0234, 0.0312, 
respectively. 

On the other hand, solutions p.6p are always unstable and are physical for: 



^fcirl^^' (3-8) 

what is equivalent to the condition T\\ < 1. 

3.2 Birth of stable 2 - cycles and transition to chaos 

Let us note that 2 - cycles are created from fixed points such that Ti e [0, 1] , Ti+i G 
[k,k + 1] with k = 1, 2, 3, 4. Therefore equations defining 2 - cycles read: 

7r (f,+i) = jY (T,) - (r,;+i - T,f + (T,+i - T,) Vr 
V,+i = -RV, + 2R (T,+i -T,)+-f{l + R)Y (f,+i^ 
Tj+i = T.i+1 - k 

^Y {f,+2) = lY (f,+i) - (r,+2 - f,+^) ' + (r,+2 - f,+i) v,+i (3.9) 
T^,+2 = -RV,+i + 2i? (t,+2 - f:,+i) + 7 (1 + i?) y (f:,+2) 

7'i+2 = T^+2 - k 
Ti+2 — Ti 

where k — 1,2,3,4 and F, Y are given by (12. 2p . We know that 2 - cycles appear 
foi" 7 = 54(TW (^^(^' " ^) + V9fc'(^'-l)' + 12(i?2 + l)2), cf. Eg.dO). 
We have solved the set of equations (13. 9p in closed form but we skip these 
lengthy formulae because of lack of space. 

The bifurcation diagram shown in Fig. 1 suggests that in the case of two fixed 
points with /s = 3,4 transition to chaos occurs after the first period doubling 



when r,i = 0. Therefore, in order to determine values of parameters at which 
the transition to chaos occurs, we have to solve Eg. p.Qp with condition Ti = 0. 
After this substitution equations p.9p are easily solved to yield: 



a^X^ + asX^ + a2X'^ + aiX + ao = 

ai^{R+lf 

a3^-k{R+l){7R+l) 

a2^k (lOfc + 3) i?2 + (lOi? - 8) P - 3k 

ai = Ak^ ((3 - i?) fc + 3) 

ao = -2fc2 (fc + 1) {2k + 1) (1 + i?2) 

Ak) _ y^X -(R+lfX+4Rk 

''='" 18(1+7?) {2+'kR.)X'-i-3iR+l)(,2k+l)X'^+6k(k+l)X+k{k+l){2k+l){R-l) 

(3.10) 
where X = Ti+i. 

For example, for R = 0.85 and fc = 4 we get from ((XTU| jir' = 0.03284 

(3) 

while for R = 0.85 and fc = 3 we obtain jcr — 0.03806 and indeed, precisely at 
this point on the 7 axis, branches of the corresponding 2 - cycles reach values 
r*i = and transform into chaotic bands. This scenario does not apply for 
R = 0.85 and fc = 1, 2 - in the case k — 2 there is another period doubling prior 
to 7cr while for fc == 1 physical solutions of (I3.10p do not exist. 

3.3 A^ impacts in one period of limiter's motion and chat- 
tering 

In the bouncing ball dynamics chattering and chaotic dynamics arise typically, 
see [17l [18] where chattering mechanism was studied numerically for sinusoidal 
motion of the table. Due to simplicity of our model analytical computations are 
possible. 

Let us assume that N = 2 impacts, Ti, Ti+i, occur in the same period. Then 
the solution A^+i = of equation (|2.1ap is always present. We thus obtain from 
Eqs. ((2ia and (l2:2a|) : 

A.+i = 0, ^(3G-2-6GT,±/d) , (3.11a) 

D = -12G'^Tl + {24:G + 12G^) T, + G^ - 12G + 16GV, + 4{3.nh) 

where G = 12-^/37, A^+i = T^+i —Ti. We can now rewrite Eqs. ()2.1|) in simplified 
form: 

Ti+i =Tj + Aj+i , , 

v,+i = -Rv, + 2i?A,+i + 3i?GA,+i (A,+i + 2T, - 1) ^^- > 



A,+i = 4^ (3G-2-6GTi + V(3G-2-6GrOVl6Gw, ) (3.12b) 

where Vi — Vi — G {iTf — 3Ti + i) is a relative velocity and the solution A^+i > 
0, cf. Eqs. (J3.11I) . was chosen. Equations p.l2p define an explicit nonlinear map. 



This formalism makes possible analysis of chattering and grazing in the TV — > cxi 
limit. It follows immediately that grazing manifold, t;* = 0, Ti+i = Ti = T*, 
i.e. Ai_)_i = 0, exists only for 3G — 2 — 6GT* < and hence for 

l>n> max {T,r, 0) , T,r = ^-^. (3.13) 

After computing eigenvalues Ai_2 of the stability matrix S on the grazing 
manifold, ?;* = 0, 3G — 2 — QGTj^ < 0, we get after straightforward calculations 
Ai = 1, A2 = i? < 1. It thus follows that the grazing manifold is always 
attracting in one eigendirection and neutral in another (T). 

In the final stage of grazing Vi is very small. Let 3G—2 — 6GTi < 0. It follows 

that ioT < Xi < I, Xi — ,_ j,^^*^"' ,ij , we get the convergent expansion of 
the square-root in ()3.12bp : 

A - -3G+2+6GTi n„ 1 ^2 \ r^ 2u 



7X1 Q X, 



+ ■ ■ ■) — -3G+2+6GTi ('^■l^j 



4G 12-^' 8-^i ' ■■•/ ~ -3G+2+6GTi 



and the inequality Xi < 1 implies: 



T,>TA,^a~^) + hl7- (3-15) 



Approximate equations describing chattering, obtained from the approxi- 
mate equation p.l4p and exact equations p.l2ap , are of form: 

ii+l = li + _3G+2+6GTi ' (3.16a) 

v,+i = Kv,, A, ^ i? (1 + p^^ji^^jj^) , (3.16b) 

and to ensure convergence ii^ — )■ we assume that Ai < 1 (i.e. Xi < |(i?^^ ~ 1)) 
- this leads to the following condition for Tf. 



T, > Ta2 = (5^3^) + y 3{l"-R)G- (3.17) 

It is now possible to estimate time of A^-th impact: 



T, 



{N) 



df ^ -r-^i+JV 






■'j=i+ 

The time T^jv) can be computed exactly since 

A _ 2tij _ 2t)j-iAj-i _ 2t'j-i r/ _ r> A /o l q^ 

"j + 1 -3G+2+6GT,- -3G+2+6G(Tj_i+Aj) -3G+2+6GTj_i ^ -"-'"j ' v^J'-^-"; 

where equations p.l4p . (j3.16bp were used, and thus 



T{N) —Ti+ _3G_|_2+6GTi Z^,^i ^ — ^i + -3G-l-2+6GTi l-fl ' (3.20) 



If the ball impacts at time Tcr < Ti < 1, performs infinite number of impacts 
and sticks (in one interval of limiter's motion) then the time is T/oo)- Solving 



the inequality ^(oo) < 1 we get conditions for grazing: 



T, > TB^^{{-2 + QG)^^{'iG + 2f~^^\>T,r, (3.21) 

T, < Tc=j^ ({-2 + 9G) + ^{3G + 2f - ^) < 1 . (3.22) 

Finally, the (approximate) condition for grazing is: max (TaiTa2,Tb) < 

T, < Tc, cf. Eqs. dSHD, (EIZI), dsnn, ^M- 
4 Summary 

It has been shown that some chaotic bands are created from fixed points after 
first period doubling in a corner-type bifurcation and critical values of control 
parameter 7 have been determined. Equations for N impacts in one period 
of limiter's motion were found and simplified significantly, making analysis of 
chattering and grazing possible. Approximate equations describing final stage of 
chattering were obtained and (approximate) condition for grazing was computed 
in analytical form. 
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